arXiv:1506.02233vl [physics.optics] 7 Jun 2015 


Propagation of time-truncated Airy-type pulses in media with 
quadratic and cubic dispersion 

Jose Angel Borda Hernandez^, Michel Zamboni-Rached^* Amr ShaarawP and loannis M. Besieris^ 

June 9, 2015 


^FEEC, State University of Campinas, Campinas, SP, Brazil 

^Department of Electrical and Computer Engineering at the University of Toronto, Toronto, ON, 
Canada 

^Department of Physics, The American University of Cairo, PO Box 74, New Cairo, Egypt 

^The Bradley Department of Electrical and Computer Engineering, Virginia Polytechnic Institute and 

State University, Blacksburg VA 24060, USA 


Abstract 

In this paper, we describe analytically the propagation of Airy-type pulses truncated by a finite- 
time aperture when second and third order dispersion effects are considered. The mathematical 
method presented here, based on the superposition of exponentially truncated Airy pulses, is very 
effective, allowing us to avoid the use of time-consuming numerical simulations. We analyze the 
behavior of the time truncated Ideal-Airy pulse and also the interesting case of a time truncated Airy 
pulse with a “defect” in its initial profile, which reveals the self-healing property of this kind of pulse 
solution. 


1 Introduction 

For a few years a new kind of localized wave solution has attracted the attention of researchers. It 
is the non-dispersive Airy pulse [T]. Contrary to other types of non-dispersive pulses, the Airy one does 
not require a spectral space-time coupling to become resistant to the dispersive effects. Instead, it just 
uses a cubic phase frequency spectrum, which greatly facilitates its experimental generation. 

An ideal Airy pulse is immune to dispersion effects of second order, but it possesses infinite energy 
making its experimental generation impossible. However, it is possible to obtain a finite-energy version 
of the ideal Airy pulse by modulating it with a time exponential function on the initial plane z = 0 pu]. 
The analysis of a finite-energy Airy pulse considering the effects of second and third order dispersion was 
first made in [3] and subsequently in [IHS]- Additional work on finite-energy Airy beams and pulses is 
provided in Refs. [Mi- 

All the previous analyses about finite-energy Airy pulses have considered exponential or Gaussian 
(initial) time apodization, but until now no one has studied the case of Airy-type pulses with truncation 
in time, that is, a temporal apodization made with a rectangular function. Such truncation, besides 
providing finite energy solutions, enables the resulting pulse to present, up to a certain distance, the very 
same characteristics of the non-truncated version. Other types of apodizations, such as the exponential 
and Gaussian, are not so efficient in this sense. 
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In this paper, we present an analytical description of Airy-type pulses truncated in time when second 
and third order dispersion effects are considered. The mathematical method used here is an extension of 
one used for describing spatially truncated Airy beams |12j . 

We analyze the behavior of the time truncated Ideal-Airy pulse and also the interesting case of a 
time truncated Airy pulse with a “defect” in its initial profile, which reveals the self-healing property of 
this kind of pulse solution. 

2 Theoretical description of Airy-type pulses truncated in time 

When quadratic and cubic dispersion effects are predominant in the evolution of a pulse in a linear 
and homogeneous material medium, the differential equation describing the pulse envelope propagation 
can be written as m 

,5C/ P2d^U ihd^U 

^ dz ~ 2 dT^ ^ Q dT^ ^ ^ 

where T is the retarded time in a frame reference moving with the group velocity {T = t — zjvg), and /I 2 
and represent the second and third-order dispersion parameters, respectively, both depending on the 
central frequency of the pulse. 

Let us define b = —^ 3 / ( 2 / 32 Tb) and the new normalized variables s = T/Tq and Z = z/ 32 /T§, where 
To is a constant that will be related with the time width of the initial pulse. 

Now, by considering an initial temporal pulse profile given by an Airy function apodized by an 
exponential, i.e, 


[/„(s, Z = 0) = Ai (s) exp (a„s) , 

with i?e(a„) > 0 , we can obtain, by solving eq.([T]), the finite energy Airy pulse given by 
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which was first obtained and analyzed in detail by Besieris et al. in [5]. Briefly, when a„ = 0 (infinite 
energy), = 0 and /32 ^ 0, we have the Ideal-Airy pulse, which is immune to the dispersion effects. In 
the case a„ = 0, /Sa 7 ^ 0 and P 2 ^ 0, the Ideal-Airy pulse can present resistance to the dispersion effects 
for long (finite) distances. 

Now, when iie(a„) > 0 the corresponding pulse given in eq.([3]) possesses finite energy content and it is 
dispersion resistant for long distances when iie(a„) >> I/Tq. However, even in these cases, the entire 
pulse starts to be destorted from the beginning and this occurs because with the exponential apodization 
at z = 0 , all time sidelobes are affected (damped) and, as it is well known, the latter are responsible for 
the pulse reconstruction. 














Based on these observations, we can envisage that a time apodization at z = 0 made by a rectangular 
function (i.e., a double step function) will preserve, within it, exactly the same initial field structure of 
the non-apodized pulse. As a result, we can expect that the resulting pulse will maintain the very same 
characteristics of its non-truncated version if the time width of the truncation is much greater than Tq. 
Actually, in the self-healing process, the inner lobes (i.e., those closer to the main peak) are fed by the 
outer ones. The latter are consummed until only the main peak remains, which then suffers distortions 
due to the dispersion effects. At this point, the pulse will have reached its maximum distance of resistance 
from the dispersion effects (depth of field). 

Next, we are going to present a mathematical method capable of describing Airy-type pulses apodized 
by a time rectangular function i.e., truncated in time, propagating in a linear medium characterized by 
both quadratic and cubic dispersion. This approach is a variant/extension of one developed to deal 
with spatially truncated Airy beams m, where the equivalent to the third order term in eq. © was not 
considered. 

Mathematically, we wish to solve (approximately) eq.(IT]), when the following initial pulse profile is 
considered at z = 0 


F{s) = Ai {s) m {s) [u {s + S) — u{s — S)] (4) 

Here, m{s) is an arbitrary function modulating the Airy function and the time truncation is repre¬ 
sented by the difference between the Heaviside unit step functions u{s + S) and u{s — S). 

We start by considering as a solution to eq. © a superposition of pulses given by 


U{s,Z) = Y. BnUn{s,Z) , (5) 

n— — co 

where the pulses C/„ are given by ©, with Bn and a„ complex constants yet unknown. 

The pulse solution ([HI) at z = 0 is written as 

CX3 OO 

U(s,0)= Y BnUn{s,0)= Y BnAi {s) exp {ans) (6) 

n— — oo n— — oo 

Our problem is to determine the values of Bn and an so that our proposed solution at z = 0 is (approxi¬ 
mately) equal to the desired truncated Airy-type pattern given by eq.(|l|). Once this is done, the resulting 
pulse propagating in the second and third-order dispersion medium is given by eq.lO. 

Let us make the following choice: 
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where a^j and L are positive constants. By using © in we get 


C/(s,0) = Ai (s) exp (aj^s) 'Y^ i?„exp(i—nsj 
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The Fourier series appearing in eq.© is suggestive. In the case we define 
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with S < -Zj/ 2, and the coefficients Bn as 


Bn = 


-s 


m (s) exp (—Ofls) exp ( —i—ns ) ds , 


( 10 ) 


the series ([9]) will represent, within the domain —L/2 < s < i/2, the following function 


A(s) = 


m (s) exp (—ors) for —S < s < S 
0 for S' < |s| < LI2 


( 11 ) 


In this way, by using (I8I9I10I11I) and with appropriate values for or and L, we can get the following 
result: 


r Ai (s)m(s) for |s| < S 

U{s,Z = Q)=l 0 for S<|s|<L/2 ^F{s) (12) 

y Ai (s) exp (ors) A (s) ~ 0 for |s| > L /2 

The suitable values for L and or serve to guarantee the result given above when |s| > L/2. Actually, as 
A(s) is the Fourier series with period L representing the function given in eg. dill) and since L/2 > S, for 
appropriate choices of L and or we have that Ai{s) exp(oR s)A(s) ps 0 for |s| > L/2 due to the behavior 
of the functions Ai(s) and exp(aRs) for Msitive and negative values of s, respectively. A good criterion 
is to choose values of or and L such thaij exp(oRL/2) >> Max[A(s)], with A(s) given in Eq. (fTTll . 

So, we have achieved our goal: The solution describing the propagation of an Airy-type pulse trun¬ 
cated in time according to eq. (na) in the presence of quadratic and cubic dispersion is given by eq. dSD, 
with On and Bn given by eqs. (TtTToI) . Or and L being appropriately chosen according the above criterion. 

We can also calculate the initial temporal spectrum of the time-truncated Airy-type pulse. Fourier 
transformation of ([ 6 ]) yields 

- ial)^ . (13) 

Examples 

Here, we are going to apply the method considering two situations where the material medium is 
fused silica, the central wavelength Aq = 1550nm (so /I 2 = —27.909ps^/km and = 0.151ps^/km), 
To = lOps, OR = 0.1, S' = 40 and L = 3S. Of course a finite number of 2N -|- 1 terms must be used in 
the summation of our solution given by eq. ®, corresponding to —N < n < N. Here we use N = 80. 

1. Ideal Airy pulse truncated in time 

Let us consider the case of an Ideal-Airy pulse with initial main peak of To = lOps, truncated by a 
finite time aperture of width forty times greater than To, i.e., S = 40. In this case we have m(s) = I and 
the pulse at z = 0 is given by eq. (HU), where the coefficients a„ and Bn are defined by the equations © 
and (Uni, respectively. As already has been mentioned, many sets of values for or and L can yield good 
results and here we have chosen L = 3S and or = 0.1. The resulting pulse is given by eq. 

Figure [1^ shows the pulse intensity at z = 0 obtained from equation ([B]), and Figure [T)d shows its 
frequency spectrum given by eq. m- 


U (0, n) = Bn exp (—exp 


Of course there are many sets of values of and L that satisfy this criterion, yielding a good result. 
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Figure 1: (a)Field intensity, at z = 0, of the Ideal-Airy pulse truncated by a finite-time aperture according 
to example 1; (b) Frequency spectrum of the truncated Ideal Airy pulse 



Figure 2: (a)Temporal evolution of the time-truncated Ideal Airy pulse intensity of example 1 at different 
propagation distances; (b) The same pulse evolution through an orthogonal projection on the plane 
{T/To,z). 
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Figure 3: Peak intensity evolution in range when S = T/Tq = 40 and 150 respectively. 


The temporal evolution of the resulting pulse at different distances is shown in Figure [5^. We can see 
that with this kind of time apodization the inner lobes, in particular the main one, are fed by the outer 
ones and maintain their amplitude and shape similarly to those of the Ideal-Airy pulse. This process 
continues until the distance (field depth) where only the main peak remains, which then is distorted by 
dispersion effects. We should note that at that distance the pulse has a shape that is very different from 
the initial one, given that there are no longer any sidelobes left, but just the main peak with the same 
time width Tq. 

For comparison, a Gaussian pulse with the same initial time and central wavelength, would have a 
dispersion length Ljj = Tq/\I32\ ~ 3.6km, while the field depth of the present truncated Ideal Airy pulse 
is approximately 35km. 

Figure shows the same pulse evolution through an orthogonal projection on the plane {T/Tq,z), 
where the accelerated character of the pulse is evident. 

It should be clear that the greater the time truncation width, the greater the depth of field of 
the resulting pulse. Figures [3Ji and [5 }d show the peak intensity evolution along the Z-direction when 
S = T/To = 40 anj] 150, respectively. Of course, a wider time width requires a larger amount of energy. 

Many other Airy-type pulses truncated in time can be described through the present method just by 
choosing the appropriate modulation function m{s) in eq. HU). 


2. Self reconstruction of a defective Airy pulse truncated in time 

Here, we are going to present the analytical description of the evolution of a time truncated Airy 
pulse with a “defect” in its initial profile. This situation can be easily described by our method, since the 
defect in question can be inserted by the function m[s) in the initial pulse profile, eq. (ED. For instance, 
let us consider the Truncated Ideal Airy pulse of the previous example with a missing piece of its initial 
field. This defect can be represented by choosing 


f 0 for —S 2 < s < —Si 
( 1 elsewhere 


(14) 


^In the case of 5 = 150 we use = 0.01. 
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Figure 4: (a)Field intensity, at z = 0, of a defective Airy pulse truncated by a finite-time aperture; (b) 
Frequency spectrum of the truncated defective Airy pulse 
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Figure 5: (a)Temporal evolution revealing the self-healing property of a defective Airy pulse truncated 
by a finite-time aperture; (b) The same pulse evolution through an orthogonal projection on the plane 
{T/To,z). 















































If this “gap” in the pulse profile is considerable shorter than the truncation time width, the Airy-type 
pulse is able to self-regenerate. To confirm this important property we set Si = 5 and S 2 = 9, such that 
the defect occupies 10% of the initial pulse. 

This defect can be clearly seen in Figldji, which shows the initial {z = 0) pulse intensity, obtained 
from eq.®. Figure Hb shows the correspondent frequency spectrum, given by ea. lfT^ . 

The temporal evolution of the resulting pulse at different distances is shown in Figl4^, where the 
self-healing property is revealed. The main part of the pulse (close to the main peak) is severely affected 
in the range 10km z < 25km, after which the main peak is regenerated, lasting until z ~ 35km, when 
it gets to be distorted by dispersion effects. 

Figure [SJd shows the same pulse evolution through an orthogonal projection on the plane (T/Tq^z). 
Again, the accelerated character of the pulse is evident. 

3 Conclusion 

Contrary to previous studies of finite-energy Airy beams and pulses, as well as of other types of 
localized waves |151I16| . based on exponential or Gaussian spatial or temporal apodization techniques, a 
novel efficient time-truncation method has been developed in this article and has been used to examine 
the evolution of a time-truncated ideal Airy pulse in a homogeneous linear medium characterized by both 
second and third-order temporal dispersion. A detailed theoretical study has been undertaken of the 
depth of field, and illustrations have been provided clearly showing the acceleration property of such a 
pulse. It has been established that a properly truncated Airy pulse can maintain its salient characteristic 
features up to ranges that are much larger than those of a Gaussian pulse with a comparable initial 
spectral structure. The same method has been used to study the evolution of a ’’deficient” ideal Airy 
pulse in the same medium, illustrating the accompanying regeneration or self-healing effects. 

The work of Chong et al. m on versatile light bullets limited to second-order temporal dispersion 
can be extended to accommodate both second and third-order dispersive effects. Eq. (HD is extended as 
follows: 


.d'tp(r,t) /32 


6 5T3 2/3o 


(15) 


dz 2 5T2 

Here, denotes the transverse (with respect to z) Laplacian operator and /3o is the wavenumber 
computed at the central angular frequency wq. Equation (USD allows a solution of the form 


= U{z,T)Q{r) 

with U{z,T) satisfying Eq.([T]) and Q{r) governed by the 3D parabolic equation 


(16) 


.dQ{r) 

dz 


= 0 


( 17 ) 


A time-truncated ideal Airy beam solution U{z,T) together with a finite-energy solution of the 
paraxial equation (HID, e.g., a Bessel-Gauss beam, will give rise to a light bullet t/’(r,t) according to 
Ea. (fT51) . In this case, the finite energy light bullet can be resistant to the diffraction and dispersion 


effects for a distance much greater than the field depth of the ordinary pulses in dispersive media. If an 
Airy beam solution to Eq. (1171) is chosen, a nonlinear transverse bending due to diffraction will appear in 


addition to the longitudinal acceleration. 
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